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Abstract 

A recent general analysis of light-baryon isospin splittings is updated and 
extended to charmed baryons. The measured S c and H c splittings stand 
out as being difficult to understand in terms of two-body forces alone. We 
also discuss heavy-light mesons; though the framework here is necessarily 
less general, we nevertheless obtain some predictions that are not strongly 
model-dependent. 
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I. BARYONS 



Mass splittings among baryons and mesons containing different valence quarks provide 
a good probe of hadronic structure and of the forces that determine this structure [|TJ. 
Recently, results were reported for a general class of models in which the baryon octet and 
decuplet splittings were assumed to arise solely from two-body interactions . It was found 
that the central masses of isospin multiplets could be fitted with a model error of 8.9 MeV. 
A similar parametrization of isospin splittings, also with a small model error, showed that 
the two measured A masses were the hardest to fit. 

In this paper, we shall extend the analysis to baryons containing heavy quarks (thereby 
severely testing the assumption that three-body forces are negligible), and also to mesons. 

Since reporting on the fit in Ref. |2j we have recognized that the two sets of input A values 
that had been averaged || were obtained using essentially the same set of input experimental 
data. Very similar analysis methods were also used, so the differences in the results should 
be treated as a systematic error. In addition, we have concluded, through discussions with 
colleagues ||, that there are additional systematic errors arising from differences in the 
treatment of the coupling to open channels. We estimate these theoretical uncertainties to 
be, at this time, about 0.5 MeV for the A states. Similar systematic uncertainties would 
also occur in E* and S* states. Since their widths are smaller and the errors are larger, we 
expect that this would give a negligible additional uncertainty to the masses of these states. 
With the enlarged error for the A masses, the fit to octet and decuplet masses has a \ 2 °f 
3.2 for 4 degrees of freedom, so that no model error needs to be invoked. However, almost 
all of this x 2 still arises from the A masses. The results from this modified fit are shown in 
Fig. |l]. The parameters from this fit will be used in fitting the heavy quark sector. 

Although a breakdown of the spectator approximation in the charmed baryon sector 
would not be at all surprising (since the heavy charm quark may significantly alter the 
environment in which the light quarks interact), we will in fact see that there appears to be 
little indication of such a breakdown in most of the data. 

The central masses of the charmed baryons may be expressed in terms of four pair energies 
S nc , T nc , S sc , and T sc (in addition to the non-charm pair energies introduced in Ref. 0) or 
in terms of the spin-dependent and spin- independent combinations D = (T — S)/4 and 
A = (3T + S)/4, respectively. (T indicates a spin-triplet quark pair, 5* a spin-singlet quark 
pair, and the subscripts give the flavor content.) However, there are only four data, so we 
introduce an additional assumption. In Ref. @ it was verified that x = D ns /D nn ~ |, as 
is often assumed in models. Here we also assume that x c = D sc /D nc ~ T7. To give this 
constraint some flexibility, we introduce the auxiliary quantity y c = \D SC — D nc , with an 
input value of 0.0 ± 1.6 MeV. Tables [I] and [Tl] summarize the fit of these pair energies to the 
A c , E c , S c , and fl c masses Including the correlated parameter errors from the non-charm 
sector shown in Table |TI| and the same model error of 8.9 MeV, we find an overall \ 2 of 0.37 
with one degree of freedom, thus the central masses are easily fitted with a sum of two-body 
energies. 

We may also express the central masses of the remaining (unmeasured) singly- charmed 



baryons in terms of pair energies to obtain the predictions shown in the left side of Table |H 



(Note that there are two S c states, since there is no symmetry restriction on the spin 
wavefunction of a use baryon. We denote the x p state by S c and the x A state by Ef . These 
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states actually mix slightly — and this was taken into account in our calculations — but 
the effects are negligible. (See also Ref. ||.) 

The fits and predictions for j = \ and j = | baryons are shown in Fig. 0. For the 
fits, the quantities plotted are the contributions of the spin-dependent terms D. The spin- 
independent terms A have been subtracted from the experimental data. For the non-charm 
baryons, we assume that D ss = D\ s jD nn . 

We have also considered an additional stronger assumption, that the A's are linear in the 
quark masses. In the non-charm sector, the three A's are then fitted with two parameters, 
with x 2 — 1-86, which is not excessive. In the charm sector, however, fitting the two new 
v4's with one new parameter gives \ 2 = 24.2. This indicates that the presence of the charm 
quark does affect the baryon structure, even though the fit to central masses shows no clear 
evidence for the existence of three-body effects. 

The fit to isospin splittings proceeds similarly. There is no advantage in discussing the 
individual E c and S c charge states, because their experimental errors are correlated in a 
complicated way. We therefore give results for the isospin tensorial quantities 

£i = £++-£°, (1.1a) 
£ 2 = £++ + £°-2£+, (1.1b) 

(1.1c) 



J c l—i c l—i c 



The new pair terms, D\ and A\ are shown in Table [TV| Note that D\ is very small and 
has little influence on the splittings; rather than treat it as a free parameter, we constrained 
it in the fit to (0.03 ± 0.05) MeV, an estimate obtained from the sizes and uncertainties of 
various terms in the simple model described in Ref. ||. The best-fit value was D\ = 0.027 
MeV. The total \ 2 f° r the fit is 6.23. The final two rows of the table are clearly the most 
intriguing: We obtain a poor fit to and a worse fit to £ 2 . If these measurements are 
correct, it appears that three-body effects are important for isospin splittings but not for 
the isospin-averaged central masses. 

The errors in the non-charm pair terms are correlated; our fit uses the full correlation 
matrix, and the shifts in the values of these terms are shown in Table |V[ The fit to the 
individual S c and S c masses gives slightly adjusted values for the central E c and S c masses, 
which were used as the input in Table |. 

As before, the pair term analysis gives rise to some predictions. These are shown in the 
right side of Table JTT[ It is interesting that the combination 2S£ L — = (—1.88 ± 0.51) 
MeV is given by — T 1 + 2T^, i.e., it is independent of the charm pair terms. Note also that 
in our model, S* 2 = £ 2 = £ 2 = £* 2 . 



II. MESONS 

A description of mesonic splittings in terms of pair energies is not useful, since a new 
parameter must be introduced for each meson. While we are thus forced to commit ourselves 
to stronger assumptions here, we still strive for results which are as model-independent as 
possible. We focus on the differential vector-pseudoscalar isospin splittings, which we denote 
by A(M). In particular, we define 
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A(B) = (B*~ - B*°) - (B- - B°) (2.1a) 
A(D) = (D*° — D* + ) — (D° — D + ) (2.1b) 
A(K) = (K*- -K*°)-(K- -K°). (2.1c) 

These combinations, which have been considered previously by a number of authors [0], are 
interesting for several reasons: 

i) A(D) has recently been measured very precisely: A(D) = 1.48 ±0.1 MeV (Ref. ||). 
This datum constitutes a severe test (or useful calibration point) for model calculations. 

n) A(K) is also known rather precisely: A(K) = -0.49 ± 0.37 MeV (Ref. §). The- 
orists have generally found it difficult to accomodate this value; even those models which 
successfully account for related splittings typically predict A(K) ps +3 MeV (see Ref. [TO]). 
It has been suggested that relativistic and/or unitarity effects (from K* — > Kit — > K*, for 
example) may be a source of the discrepancy. 

Hi) B~ - B° is known to be -0.12 ± 0.58 MeV (Ref. 0), but the B* isospin splitting 
has not yet been measured, thus there is room for a prediction of A(B). 

iv) As we shall presently see, it is possible to obtain predictions for the ratios A (B) : 
A(D) : A(K) which are quite model-independent. 

We write the quark- mass and electromagnetic contributions to A(M) as A m (M) and 
A 7 (M), respectively. To obtain A m (M), we employ an interpolation technique first sug- 
gested by Chan Consider, for the sake of definiteness, A m (D). From equation fl2.1b|) , 



we see that A m (D) m (m u — m d )(dE hyp /dm q )\ rnq=mu . This latter quantity may be accu- 
rately estimated with very little reliance on theory, since the function Eh yp (m q ) is rather 
well-determined by measurements of D* — D, D* — D s and J/^f — r] c ; an interpolating curve 
drawn through these points allows the extraction of (dEh yp / dm q )\ mq=rriu [Q. The same 
procedure can be applied to A m (B) and A m (K), and the relevant curves for all three cases 
are shown in Fig. [|. 

In the case of D and B mesons, the hyperfine splittings turn out to be surprisingly 
independent of the lighter quark's mass: 

D* - D = 141.4 ±0.1 MeV, 
D* s -D s = 141.5 ± 1.9 MeV, 
B* - B = 46.0 ± 0.6 MeV, 

B* - B s = 47.0 ± 2.6 MeV . (2.2) 

(This is qualitatively understood in the quark model, where hyperfine splittings vary with 
the quark masses as |-?/'(0)| 2 /m 1 m2; in a linear potential, (^(O)) 2 is proportional to the 
reduced mass fi, so that Eh yp oc (mi + m2)~ 1 , which becomes independent of mi in the limit 
772,2 >> mi.) Thus A m (D) and A m (B) are both approximately zero. The second column of 
Table [VT] shows the actual results of the interpolation procedure; we performed a quadratic 
fit to the three data points for D mesons, and a linear fit to the B* — B and B* — B s data 
points for B mesons. 

The strange meson sector differs in that Eh yp has significant m^-dependence: the available 
data here are K* — K = 398 MeV, D* - D s = 142 MeV, and B* - B s = 47 MeV. A cubic 
spline interpolation gives (dEh yp /dm q )\ mq=mu = —0.36, however, without a data point for 
the ss mesons, the distance over which we must interpolate in this case is unsettlingly large. 
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The ss data point is not directly available, of course, because annihilation forces rotate the 
pseudoscalar ss meson (r] s ) and its nonstrange partner (rj n ) into the physical 77 and 7/ states. 
Nevertheless, we may estimate the mass that r] s would have in the absence of such forces by 
making use of the empirical rule |13] for the difference of like-flavor vector and pseudoscalar 

0.56 GeV 2 . 



squared masses: m v — m P m 0.56 GeV^. This implies m na = \Jrri^ — 0.56 = 690 MeV. 
Independently, an equal-spacing rule should hold for the squared masses of ideally-mixed 
pseudoscalars so that m 2 a = 2m 2 K + m 2 , which gives m Vs = 687 MeV. Using this additional 
data point gives d(Eh yp / dm q )\ mq=mu = —0.32 ± 0.06, where the error estimate reflects the 
dependence of the fitted slope on our choice of constituent quark masses [Q, and on the 
value of m Vs . 

Almost all model calculations are consistent with — m u = 5 ± 1 MeV We adopt 
this value and show A m (M) in the third column of Table |VT[ 



The electromagnetic contribution to A(M) is given by, 

A W = a . m ,J/i\ + £!*!^-> , (2.3) 




r I q 

where the subscript h refers to the heavy-flavor quark (h = s,c,b), and the subscript n 
refers to the light anti-quark. The quantity a e m<?/(l/ r )o (1) is the Coulomb energy in the 
pseudoscalar (vector). The difference (1/V) — is induced by hyperfine distortion 

of the wavefunctions; though it is often assumed to be negligible, explicit potential-model 
calculations indicate that this Coulomb term is only about 30% smaller than the magnetic 
term in (|2.3j ) (We find it most convenient to discuss A 7 (M) in a potential- model 



framework, however, our results will turn out to be nearly independent of the parameters 
- and even the functional form — of the potential, hence we suspect that our conclusions 
are more reliable than our method.) 

Consider calculating A 7 (M) in a Hamiltonian whose spin-independent part is 

H st = ^- + br-^(-\ . (2.4) 



2fi 

This is a standard Coulomb-plus-linear Hamiltonian, with a particular ansatz for the way 
in which the effective strong coupling constant varies with the reduced mass of the system: 
a s {ii) = ot° s ■ (b/fi 2 ) 1 / 3 . We shall presently show that this form for a s (fi) gives excellent 
phenomenology, but for the moment we just capitalize on the fact that it yields simple 
scaling laws for the wavefunctions and energies of H si . Since 

/i^ = ^ + (^)r-^(M) V3 , (2.5) 

it follows that the energies are proportional to (6 2 //x) 1//3 , and wavefunction properties with 
dimension (length)" are proportional to (fib)~ n / 3 . Upon adding the spin-dependent term, 

H hyp = 3 -^a s ^)^^S 3 (r), (2.6) 
9 m h m n 
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we find the following scaling laws for the hyperfine, magnetic, and Coulomb energies: 



32tt . JV(0) , v fib 

E hyp = __ a ,(/i)^-^- oc a s (n) — — (2.7) 
9 m h m n m h m n 

87r|^(0)| 2 ujb , < 

Emag = OL em q h — CX (2.8) 

3 m h m n m h m n 

„ __ r/l\ 64tt oJjm) ^ (0|<5 3 (r)|n)«|0) fib 

Ecoui = a em q h U-) -{-) \ = — 2^ B ET^ — « ■ 2 - 9 

Thus the Coulomb and magnetic terms scale in the same way, and we have 

&y{B) : A 7 (L>) : A 7 (K) = q b Eh ^ B ) : q ^^} : q* ^*} . (2.10) 

a>(/ifl) OsK) 

In fact, it is easy to check that the result (|2.10| ) holds not just for the specific potential 
in (|2.4j) , but for any potential of the form 

w = Cr ,_^(JL) . (2 . n) 

(And also for V(r) = Clog(r/r ) — (4/3)(a°/r)(C//i) 1 / 2 , which behaves in many ways like 
the (3 — > limit of Vp{r)\ we will refer to this logarithmic potential as Vo (/").) In the general 
case, the ratios of electromagnetic splittings are 

A 7 (5) : A 7 (L>) : A 7 (#) = q b fi e B E hyp (B) : q c ^ D E hyp (D) : q s ^ K E hyp (K) , (2.12) 

where e = (/? + + 2). 

With regard to the phenomenological soundness of equation ( |2.11| ), Fig. f| shows the 
"running" of ot s (fj) for three cases of interest, namely (3 = 0, 1, 2 (corresponding to logarith- 
mic, linear, and harmonic confinement, respectively). In each case, the curve lies close to 
the empirical values of a s (fj, K ), a s (fJ, D ), and a s (fj, B ), which were determined by fitting to the 
measured 2P — IS and 3 Si — 1 S meson splittings [jnj. (Of course, the success of these fits 
is not mysterious; there is little difficulty in accurately parametrizing a monotonic, closely 
spaced set of data points (fx g —/jl k is only about 0.1 GeV) with an adjustable fitting function 
which is similarly monotonic.) Thus the potentials Vg(r), originally constructed just to give 
simple scaling behaviour, actually coincide closely with purely phenomenological potentials 
tuned to fit the meson spectrum. This indicates that we may apply equation ( [2.12j ) with 
some confidence: fitting A 7 (L>) to A(Z>) expt - A m (D) = 1.50 ± 0.12 MeV, and taking (3 = 1 
we obtain the results shown in columns 4 and 5 of Table |VI[ The results are quite insensitive 
to (3 and to the constituent quark masses; A 7 (B) and A 7 (.fr) change by less than 5% when 
beta is varied between and 2, and by less than 10% when the quark masses are varied flTJ 



We thus assign a 10% "parameter uncertainty error" to each of A 7 (B) and A 7 (i^), and add 
it in quadrature with the error arising from A 7 (.D). 
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Both of the predictions in Table [VI] are interesting. First, we obtain a negative value for 
A (5), in disagreement with most other authors (see Table |VT|). The sign of our prediction 
follows just from the near- vanishing of A m (B) (as implied by the near-equality of the B and 
B s hyperfme splittings), and the fact that A 1 {B) is negative definite (since the Coulomb 
and magnetic terms in equation ( |2.3| ) have the same sign). Second, our result for A(K) is in 
agreement with the measured value, and even has the correct sign (though in this case there 
is a large cancellation between the quark-mass and electromagnetic terms), so that A(K) 
can apparently be explained without invoking large relativistic or unitarity effects. 



III. CONCLUSIONS 



In summary, we have analyzed baryon and meson isospin splittings and isomultiplet 
central masses in a constituent quark picture, avoiding (as much as possible) strong reliance 
on particular models or parameter sets. Our results indicate that improved measurements of 
the A, E c , H c , and B* splittings would be most interesting. The A masses are currently the 
most discrepant in the non-charm sector, the S c and the S c seem to indicate the existence 
of substantial three-body effects, and measurement of the sign of A(£) = (B*~ — B*°) — 
(B~ — B°) will provide a significant test of model calculations. 
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FIGURES 



FIG. 1. Octet and decuplet isospin splittings (deviations from multiplet centers) in MeV, for 
multiplets with / > 0. The fitted masses of j = \ states are indicated by filled circles, and the 
fitted masses of j = | states are indicated by open circles. The crosses give the experimental data 
and errors. For clarity, 3 = \ and j = | points have been slightly offset. 

FIG. 2. Multiplet hyperfine energies, in MeV. A spin-independent mass term that depends on 
the constituent quark flavors has been subtracted from each energy. The fitted masses of j = ^ 
states are indicated by filled circles, and the fitted masses of j = § states are indicated by open 
circles. The data and errors are shown by crosses, but only the S c and f2 c have errors of noticeable 
size. 

FIG. 3. Hyperfine splittings in sq, cq, and bq mesons, as a function of the antiquark mass. (See 
the text for an explanation of the ss data point.) Typical interpolating curves are also shown. 

FIG. 4. The data points show effective values of a s (fJ> K ), a s (fi D ), a s (fi B ), obtained by fitting 
meson splittings to harmonic+Coulomb, linear+Coulomb, and logarithmic +Coulomb potentials. 
The curves show the corresponding behaviour of a s (fi) predicted by equation ( |2.11| ). 
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TABLES 



TABLE I. Contributions to the fit to charmed baryon central masses. 





T 


T 

± ss 


Snn 




A 

■"-■tic 




D nc 


D sc 


total 


data 


± 


x 2 


A c 








213.5 





2068.7 











2282.3 


2284.9 


0.6 


0.053 


S c - A c 


412.1 





-213.5 











-27.0 





171.5 


167.8 


0.3 


0.047 













350.6 


1034.4 


1088.5 








2473.3 


2470.0 


1.5 


0.130 


n c 





557.0 











2176.9 





-17.6 


2716.3 


2719.0 


7.4 


0.054 


Vc 




















-6.8 


6.6 


-0.15 


0.00 


1.60 


0.008 



TABLE II. Comparison of the non-charm pair terms obtained from fits to charmed and 
non-charmed baryon central masses. 





T 


T 


T 

- 1 SS 




Sns 


From non-charm data 


412.2 ± 2.6 


486.8 ± 3.0 


556.8 ± 2.8 


212.5 ± 5.4 


351.4 ± 4.4 


Including charm data 


412.1 


486.8 


557.0 


213.5 


350.6 


x 2 


0.002 


0.000 


0.005 


0.034 


0.033 



TABLE III. Predictions for charmed baryon masses and isospin splittings. 

Baryon Predicted mass (MeV) Splitting Prediction (MeV) 

S* 2494 ± 16 1.00 ± 0.52 

E' c 2587 ± 12 S* 2 1.64 ± 0.21 

E* 2621 ± 15 E'} -0.52 ± 0.33 

n* 2743 ± 17 Ef -0.44 ± 0.37 



TABLE IV. Contributions to the fit to charmed baryon isospin splittings. 





T 1 


rp2 


si 


Al 


Dl 


total 


data 


± 


x 2 




-1.29 








2.24 


-0.11 


0.84 


0.90 


0.42 


0.018 







1.64 











1.64 


-1.30 


1.48 


3.95 










-4.03 


1.12 





-2.91 


-6.30 


2.30 


2.17 
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TABLE V. Comparison of the non-charm pair terms obtained from fits to charm and non-charm 
isospin splittings. 

T[ Si 

From non-charm data -1.29 ± 0.0 1.70 ± 0.21 -1.60 ± 0.26 -4.02 ± 0.10 

Including charm data —1.29 1.64 —1.59 —4.03 

X 2 0.000 0.082 0.000 0.006 



TABLE VI. Quark mass and electromagnetic contributions to A(D), A(_B), and A(K). 



Meson 


[dEfoyp / dfTlq) \ mq=mn 


A m (M) 


AT(Af) 


A(M) 


D 


0.004 ± 0.01 


-0.02 ± 0.06 


1.50 ± 0.12 


1.48 ± 0.10 a 


B 


0.005 ± 0.01 


-0.02 ± 0.06 


-0.26 ± 0.03 


-0.29 ± 0.07 


K 


-0.32 ± 0.06 


1.60 ± 0.44 


-1.75 ± 0.22 


-0.15 ± 0.49 



a Used in fit. 



TABLE VII. Various predictions for A(B). (The fourth entry is actually a constraint obtained 
by fitting the measured decay constants /d and fs to a nonrelativistic quark model.) 



Source 
Ref. |lP 
Ref. U| 
Ref. |9| 
Ref. |0| 
This work 



A(g) 
0.3 ± 0.03 
0.1 
0.3 

-0.05 to 0.49 
-0.29 ± 0.07 
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